Spin-accumulation and Andreev-reflection in a mesoscopic ferromagnetic wire 
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The electron transport though ferromagnetic metal — superconducting hybrid devices is con- 
sidered in the non-equilibrium Green's function formalism in the quasiclassical approximation. At- 
tention if focused on the limit in which the exchange splitting in the ferromagnet is much larger 
than the superconducting energy gap. Transport properties are then governed by an interplay be- 
tween spin-accumulation close to the interface and Andreev reflection at the interface. We find that 
the resistance can either be enhanced or lowered in comparison to the normal case and can have 
a non-monotonic temperature and voltage dependence. In the non-linear voltage regime electron 
heating effects may govern the transport properties, leading to qualitative different behaviour than 
in the absence of heating effects. Recent experimental results on the effect of the superconductor 
on the conductance of the ferromagnet can be understood by our results for the energy-dependent 
interface resistance together with effects of spin-accumulation without invoking long range pairing 
correlations in the ferromagnet. 
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I. INTRODUCTION 

Much theoretical and experimental work has addressed 
the effect of a superconductor (S) in proximity to a nor- 
mal metal (N) on the transport properties during the last 
years, see Ref. [j] and references therein for a review. Most 
experimental results can be explained in the framework 
of the quasiclassical theory of superconductivity account- 
ing for a "long range" proximity effect with a coherence 
length £ = (hD /2/cbT 1 ) 1 / 2 , where D is the diffusion coef- 
ficient of the normal metal and T is the temperature. On 
the other hand, applications of the quasi-classical theory 
to transport in heterostructures containing ferromagnets 
(F) are still scarce. In contrast to normal metals the 
presence of a strong exchange field in the ferromagnet 
leads to a strong difference in the energy dispersions for 
the two spin bands. However, long-range coherence in 
normal metals requires spin degenerate bands close to 
the Fermi energy, since singlet superconductivity couples 
quasiparticles of different spins by Andreev reflection. 
The consequence of the exchange field energy h xc is a 
strong decoherence of quasiparticles belonging to the dif- 
ferent spin bands. Typically the superconducting energy 
scale A is smaller than h xc by several orders of magnitude 
for (Al, Nb) vs. (Fe, Ni, and Co), respectively. Thus, 
the proximity effect in ferromagnetic metals is negligible 
and a ferromagnet in contact to a superconductor may 
be considered as an incoherent metal coupled to the su- 
perconductor. In this case all changes induced by the 
contact to a superconductor depend on the properties of 



the interface itself n |This is accomplished by the effect of 
spin accumulationntj, which requires no phase coherence 
in the ferromagnet and can therefore have a much longer 
range than the proximity effect. The main purpose of 
this paper is to study the mutual influence of resistance 
changes by spin accumulation and interface properties. 

Recently heterostructures of ferromagnets and su- 
perconductors phave been experimentally realized and 
investigated.ET0 Several unusual phenomena have been 
unveiled. riThe experimental results in point contact 
geometriesQ can be explained by the reduced, bias- 
dependent transparency of the interface due to spin- 
dependent band mismatch between the normal metal 
and the ferromagnet na The .ex-perimental results in diffu- 
sive nanostructured samplesErLl are more intriguing. The 
measured conductance changes on the ferromagnetic side 
can be positive and negative at the superconducting tran- 
sition with amplitudes much larger than anticipated. The 
sign and the amplitude of the changes appear to depend 
strongly on the ferromagnetic - superconductor interface 
transparency. It has been conjectured that a strong mu- 
tual influence of the superconductors and ferromagnetic 
conductors and a penetration of the superconducting or- 
der parameter into the ferromagnet over distances many 
times longer than expected frftift the above estimates 
might explain the observations.Erlj 

Some effects of the interplay between spin accumu- 
lation and Andreev reflection in diffuse systems have 
been discussed in Refs. ^|. Since the spin-current into 
a superconductor vanishes at sufficiently low bias and 
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temperature, a non-equilibrium spin-accumulation builds 
up on the ferromagnetic side in order to conserve the 
spin-currents. The spin-accumulation causes an addi- 
tional boundary resistance which is of the order of the 
resistance of the ferromagnetic wire of a length of the 
spin-flip diffusion length. Therefore the resistance of 
the F-S system should be always larger than that of 
the F-N system, in contradiction with many of the ex- 
perimental observations. A possible reason for this ap- 
parent failure is the neglect of changes in the interface 
resistance in the transition from F-S to F-N. Previous 
theories took into account only perfectly ballistic inter- 
faces for which the resistance is determined purely by 
the matching of the adjacent Fermi surfaces. The in- 
terface resistance and its modulation are of the order of 
the Sharvin resistance, which is negligible compared to 
the total one. However, in thej-sputtered samples with 
relatively large contact areasoQ the interface can con- 
tribute significantly, especially when differences of resis- 
tances below and above the superconducting transition 
temperatures are considered. Other transport phenom- 
ena in ferromagnet-superconductors systems have been 



studied in Ref. 10 



It has been speculated that the triplet component of 
the order parameter induced by the fluctuations of the 
spin-orbit-scattering potentials is essential in mesoscopic 
junctionsJlJ Neglecting magnetic impurities and spin- 
orbit coupling the superconducting order parameter is 
a spin singlet. However, magnetic impurities or spin- 
orbit coupling, induce a fluctuating spin triplet compo- 
nent with zero average. The triplet component is 'long- 
range' coherent in the ferromagnet since it couples elec- 
trons and holes with the same spin and the exchange 
field in the ferromagnet does not play a role. However, 
the contribution to the conductivity from the triplet fluc- 
tuations is only relevant when the fluctuations are rela- 
tively large which is only the case when the conductance 
is close to J-he quantum conductance. The experimen- 
tal samplestrQ have a much larger conductance, and we 
do not expect that such mesoscopic fluctuations play an 
important role. 

None of the above-mentioned theories can explain the 
recent experimental results. This makes it necessary to 
study the properties of the contact between the ferro- 
magnet and the normal metal in more detail and to ac- 
count for a possible spin accumulation and heating ef- 
fects in the ferromagnet including all possible interfaces 
between the ferromagnets and the superconductors. In 
particular we will go beyond tbje assumptions of a perfect 
transparent metallic interfaces and discuss its influence 
on the observed conductance changes below the super- 
conducting transition temperature and for bias voltages 
less than the superconducting gap. We will in this work 
radically disregard the proximity effect. Therefore our 
results only apply to ferromagnets with h xc 3> A, which 
is e.g. the case for the magnetic transition metals (Fe, 
Co, and Ni) in conjunction with superconducting met- 
als like Nb and Al. This assumption is supported by 



the experimental fact that FS interferometers show no 
phase-periodic oscillations down to the level of 0.1e 2 /h 
in strong ferromagnets a~u In contrast to the calculations 
presented in this paper, the proximity effect could be im- 
portant in weak ferromagnets. Below we will show that 
most of the recent experimental results can be explained 
in terms of the energy-dependence introduced by inter- 
face conductance and the accomplishing change in the 
spin accumulation. It is important to note that these 
changes are small in comparison to the total resistance, 
which is dominated by the long ferromagnetic wire. Nev- 
ertheless they play a dominant role when only the resis- 
tance changes are measured. 

The paper is organized in the following way: Section 
H gives a description of the diffusive ferromagnetic wire 
both in the limit of elastic and inelastic scattering be- 
tween the electrons. Section III treats the boundary con- 



dition between the ferromagnet and the superconductor 
which is crucial to the understanding of the transport 
properties. The results for the conductance obtained 
from the description of the ferromagnetic wire with the 
boundary conditions are discussed in Section IV . Finally 
we compare our results with experiments in Section |V] 



and give our conclusions in Section VI 
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wire. It is placed between a ferromagnetic reservoir held at 
voltage V and a grounded superconducting reservoir. The 
contact to the superconductor is an arbitrary connector, char- 
acterized by spin-dependent conductances in the normal state. 



II. DESCRIPTION OF THE FERROMAGNET 

We consider a ferromagnetic diffusive wire connected 
to an ideal (ferromagnetic or normal metal) reservoir on 
one side and to a superconducting reservoir on the other 
side as depicted in Fig. ^. The wire is characterized by 
length L, cross-section A and spin-dependent conductiv- 
ities <7f and crj,. In this Section we discuss the kinetic 
equations describing the ferromagnetic wire in the ab- 
sence of the proximity effect. We consider collision with 
impurities to be the dominant scattering process and use 
the diffusion approximation. The electrons in the quasi- 
one-dimensional wire are described by energy e and spa- 
tial x dependent distribution functions f s (e,x) for the 
two spin directions s = +, — for spin f and J., respectively. 
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The distribution functions obey two coupled Boltzmann 
equations in the diffusive limit. Other scattering mecha- 
nisms will be specified in the following subsections. 

Instead of the spin-dependent conductivities, it is con- 
venient to introduce the total conductivity a = <jf + cr^ 
and the spin-polarization of the conductivity 7 = (cj — 
<7jJ/<7. The spin dependent conductivities are then ex- 
pressed as a s = (1 + s"f)a/2. We will also make use of 
the total conductance (resistance) of the ferromagnetic 
wire G F = Aa/L (R F = l/G F ). 



A. Elastic scattering 

In the elastic scattering case the energy is conserved 
in the scattering processes. This makes it necessary to 
study the energy dependent distribution functions in the 
ferromagnetic wire. In addition to elastic impurity scat- 
tering we consider only the spin-flip scattering processes, 
accounted for by the spin-flip length l s f. Then, the kinetic 
equation reads 

' /2 f s (e,x) = ±-[f s (e,x)-f- s (e,x)] . (1) 



dx 



sf 



The current for spin s is given by (e = |e|) 
A f , df s (e,x 



A f df s (e,x) f 
Is{x) = - as - J de —— = J 



del s {e,x). (2) 



This equation defines the spectral current I s (e,x). Elec- 
trical and spin currents are /charge = I] (x) + I± (x) and 
/spirit) = I](x) — I[(x) and similar for the spectral cur- 
rents. It is convenient to introduce the conductivity- 
averaged distribution function 



(7t 1 

f e l(e,x) = — fi(e,x) + —fi(e,x) 
a a 



(3) 



and the nonequilibrium spin distribution function 

f sv {e,x) = Me,x) - fi(e,x) . (4) 
The kinetic equations in terms of these functions read 



d 2 

—^f e l(e,x) = 

d 2 1 

^2/s P (e,2;) = p"/ S p(e,x) 



(5) 
(6) 



The first equation is the spectral current conservation 
and the second describes spatial relaxation of the non- 
equilibrium spin-distribution. 



B. Inelastic scattering 

The reason to investigate the role of inelastic scatter- 
ing is the convenient fact that the ferromagnet is an 



incoherent metal with rather strong correlations. Both 
phonon and electron-electron scattering can mediate in- 
elastic scattering. In general it is not obvious which 
should dominate and both should be treated on equal 
footing. In order to achieve insight in the physics it is 
useful to consider limiting cases as well. 

In the limit of strong inelastic scattering we assume 
that the electron-electron interaction is stronger than the 
electron-phonon relaxation. When a bias voltage is ap- 
plied, the local electron temperature can therefore be dif- 
ferent from the temperature in the reservoirs. This trans- 
port regime is relevant when the typical inelastic scatter- 
ing length is smaller than the spin-diffusion length. 

The electrons relax to a local equilibrium 



fs{t,x) = f(e;ii s (x),T e i(x)) , 



(7) 



where (J> s (x) is the spin-dependent chemical potential, 
T c \(x) is the local temperature and 



1 + exp((e - n) /ksT) 



(8) 



is the Fermi-Dirac distribution function. This makes it 
possible to integrate the kinetic equation and the cur- 
rents over energy and to obtain equations for the local 
chemical potentials and electron temperature. 

The spin-dependent (electric) current from Eq. (^) is 



Current conservation requires 



— (<x T fi t (x) +a l ni(x)) = 0. 



(9) 



(10) 



Spin relaxation occurs within the spin-diffusion length 

d 2 1 

[fif(x) - m(x)] = p- [fif(x) - m(x)] . (11) 

The local spin-dependent chemical potentials in the ferro- 
magnet are determined by ( |l0| ) and (|ll]) and the bound- 
ary conditions to be discussed below. 

Additionally, we need equations describing energy 
transport in the system to account for heating of the 
electrons. The energy current is 



U*) = ~?2s°' J 

s 

= [^(x)I t (x)+ f i l (x)I i (x)}/e + lQ(x), (12) 
where the heat current is 

dT e i(x) 



Iq( x ) = -Kq(x)A- 



dx 



(13) 



the heat conductivity kq(x) = crCoT e \(x) and the Lorentz 
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number is £ = \ ( ^f- 
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The conservation law for the energy current dictates 



dt 



(14) 



where p € (x) is the local energy density. The energy re- 
laxation between the electronic system and the phonons 
at sufficiently low temperatures ialj 



dt 



C 



rel. 



{k B T) 5 - {k B T el (x)Y 



(15) 



where £ parameterizes the strength of the electron- 
phonon interaction C = 48Tr((5)N(e F )\* /(h 3 LU 2 D ), C(5) » 
1.04 is the Riemann zeta function, N(cf) is the density of 
states of both spins per unit cell, A* is of the same order 
of magnitude as the electron-phonon coupling constant 
A, and Hu)£> is the Debye energy. 

The conservation of energy current (04) together with 
the expression for the energy relaxation (|15|) give a differ- 
ential equation for the local electron temperature which 
can be solved together with the boundary conditions to 
be discussed below. 

When the electron-phonon interaction is weak there 
is no exchange of energy between the electron and the 
phonon systems so that the right hand side of Eq. Jl4| ) 
can be set to zero and we have conservation of the energy 
current due to the electron transport dl e {x)/dx = 0. In 
the opposite limit of strong electron-phonon interaction 
the electron temperature equals the lattice temperature. 
The differential equation for the energy conservation with 
the boundary conditions given above can in these two 
cases be solved exactly. In the intermediate regime the 
equations will be solved numerically. 



III. BOUNDARY CONDITIONS 

The condition that the ferromagnet should be com- 
pletely incoherent leads to simplified boundary condi- 
tions for the kinetic equations. These boundary con- 
ditions can be derived from the boundary conditions 
for the quasiclassical Green's function.^! A transparent 
form suitable, for diffusive systems has been presented 
by Nazarov.Ej We will follow the spirit and the notation 
of this paper. A circuit theory for ferromagnetic-normal 
metal systems has been presented in Ref. [15| A contact 
is described by a set of transmission eigenvalues {T n } or, 
equivalently, by a distribution of the transmission eigen- 
values p{T). The boundary condition at the contact is 
expressed through a conservation law for the matrix cur- 
rent in the Keldysh formulation. In the framework of 
superconductivity it is a 4 x 4-matrix comprising 2x2 
Keldysh space and 2x2 particle-hole (Nambu) space. 
The matrix current through the FS-contact istJ 



4l> 



(GfGs — GsGf) 



4 + T n (G F G s + G s G F -2) 



(16) 



This matrix current has to be equated to the diffusive ma- 
trix current entering the contact from either side. The 
two sides of the contact are characterized by the Keldysh 
matrix Green's functions Gs and Gf, which we will spec- 
ify to be the superconducting reservoir and the ferro- 
magnetic wire, respectively. The Keldysh-Nambu matrix 
Green's function of the superconductor in equilibrium is 



Gs{£) -{ Gg (e) 



(17) 



A similar structure holds for any matrix in Keldysh 
space. In local equilibrium the Keldysh (1,2) component 
in Nambu space is 



Gf(e) = (Gf(e)-G£(e))(l-2/ s (e)) 



(18) 



where / s (e) = [1 + exp(e/fcBT)] _1 is the quasi-particle 
distribution function in the superconductor and we have 
set the chemical potential in the superconductor to zero. 
Gg (e) and Gg-(e) are retarded and advanced Nambu 
Green's functions determining the spectral properties of 
the superconductor. In the BCS case with a real order 
parameter they are 



G R (e) 



G A (e) 



(e + i0)f 3 - zAfi 



(19) 



The diagonal component represents the normal retarded 
Green's function whereas the off-diagonal component is 
conventionally called the anomalous Green's function. 
On the ferromagnetic side we completely neglect the 
proximity effect leading to the spectral functions Gp = 
T3 = — Gp. The absence of an anomalous component 
is a result of the absence of the proximity effect. The 
Keldysh component accounts for the spin-dependent non- 
equilibrium distribution: 



1 



2/ T F (e) 






2/f(-«) 



(20) 



where /^(e) and ff(e) are the quasi-particle distribu- 
tion functions close to the interface on the ferromagnetic 
side. The spectral electrical current is determined by the 
Keldysh-component of the matrix current according to 



l el (e) = -Tr 



(21) 



Eq. (^0|) and Eq. ( pi] ) suggest a representation of the 
diagonal components of the Keldysh component of the 
matrix current in the form 



(22) 



Now we are in the position to calculate the spin-resolved 
currents through the contact. Performing the calcula- 
tions along the lines of Ref. [l4| we find the spectral spin- 
dependent current 



4 



Is(e) = 



Gqp(£) (f S (e)-ff(e)) 



+ 



2e 
G A (e) 
4e 



(23) 



(l-/f(e)-/^ s (-e)) 



The quasiparticle conductance GQp(e) and the Andreev 
conductance G^(e) are determined by the properties of 
the contact and the spectral properties of the two metals 
connected by the contact. The distribution of transmis- 
sion eigenvalues can be incorporated in a single charac- 
teristic complex function 



Z[x) = - 



nh^ 2 + T„(x-l)' 



(24) 



where x(e) = Tr{Gf (e), Gp 1 (e))}/4. The conductances 



Gqp(c) = ReZ{x)Rex- 

G A (e) = -^|l 
tmx 1 



ImZ(x) 
Imx 

„,2| 



Im 2 Vl 



(25) 
(26) 



The contact is characterized by a transmission dis- 
tribution, which leads to contact-specific energy depen- 
dences of the conductances. The normal state conduc- 
tance is Gbn = (e 2 /27r?i) J2 n Tn- For a ballistic model 
contact all transmission eigenvalues are equal to one 
for the propagating channels and zero otherwise and 
— N, where N is the number of propagating 
channels. The distribution function in the case of a dirty 
interface isli3 



p(T) = -^Gbn 



1 



T 3 / 2 v / T :r T 



(27) 



and in the case of a diffusive contact the distribution is0 

(28) 



P (T) = A Gbn _L= 



Finally, for a tunnel conductance a perturbation expan- 
sion in terms of the small transmission eigenvalues can 
be performed. We list the characteristic function Z{x) 
for a number of generic contacts in Table (|): Tunnel 
junction, ballistic contact, diffusive contact and dirty in- 
terface. In the case of an incoherent metal on one side 
(i.e. Gp = fa), the argument of the characteristic func- 
tion reduces to x — Ttt^G^ /2. The result in this case 
is demonstrated explicitely in Table | for a contact of a 
BCS-superconductor with spectral functions (|l9|). The 
energy dependence of these spectral conductances is de- 
picted in Fig |2[ Below the superconducting gap only the 
Andreev conductance is nonzero, gradually decreasing 
from the value of 2Gbn for the metallic junction to zero in 
the tunnel junction. Above the gap the Andreev conduc- 
tance vanishes rather quickly ~ 1/e 2 - Also quasiparticle 
transport becomes possible and, thus, spin-transport into 
the superconductor. 

The properties of these contacts are demonstrated by 
the temperature dependence of the linear conductance 
following from 



Gbs(T) = de(G QP (e) + G A (e)) 



df(e,0,T) 
de 



(29) 



This is the conductance that would be measured if the 
contact would be placed between a normal reservoir and a 
superconducting reservoir. The temperature dependence 
of the contact conductance (^) is shown in Fig. [|. The 
dashed and dotted lines show the conductance of the dif- 
fusive contact and the dirty interface, respectively. The 
resistance of the diffusive contact shows the well known 
reentrant behavior, i.e. it reaches the normal state con- 
ductance at zero temperaturc.li3 The resistance of the 
dirty interface after a small drop below the critical tem- 
perature is higher than the normal state value and satu- 
rates at low temperature at v^Abnm 
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abbreviations 



x = |Tr{Gs , Gp} 



TABLE I. Spectral conductances of different generic contacts defined by the transmission distributions in the second column. 
The characteristic function Z of the contact was defined in Eq. ( |24| ) . As an example we present electrical and spin conductance 
for a contact between a ferromagnetic metal and a BCS Superconductor in the last two columns. For energies below the 
superconducting gap A the quasiparticle conductance Gqp vanishes for all contacts. Note, that the energy argument of all 
quantities e is understood to be the absolute value of the energy. 
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FIG. 2. Spectral conductances for different types of con- 
tacts. The solid curves denote the quasiparticle conduc- 
tance Gqp(e) and the dashed curves the Andreev conduc- 
tance Ga(c)- The contact types as indicated in the figure are 
metallic junction (all T n — 1), diffusive contact (transmission 
distribution as defined in Eq. (53)), dirty interface (Eq. (5?])), 
and tunnel junction (all T n <C I). 



Additionally, we introduce mixed contacts as a model 
for an inhomogeneous interface with distributed regions 
with low and high transparency. The relative admix- 
ture q of a tunnel and (1 — q) of a ballistic contact al- 
lows switching continuously from one limit to the other, 
covering approximately the universal cases of a diffusive 
contact (q w 0.5) and a dirty interface (q ss l/y/2 with 
a single parameter q. A common feature of the temper- 
ature dependence of all these contacts except the tunnel 
junction is that right below T c the resistance drops. At 
lower temperatures the resistance increases again except 
in the case of the purely ballistic contact. The drop of re- 
sistance of these contacts close to T c can be traced back 
to the temperature dependence of the superconducting 
order parameter A(T). The resistance drop is caused by 
the leading order contribution of the change in the su- 
perconducting gap A(T) oc (1 — T/Tf) 1 / 2 to the Andreev 
contribution and the conductance.Ej 

The boundary condition presented so far imply that 
the transmission ensembles and the number of channels 
are the same for the two spin species. In reality the trans- 
mission matrices for spin-up and spin-down states can 
be different. A microscopic calculation of the transmis- 
sion eigenvalues is beyond the scope of the present paper. 
We will therefore heuristically generalize the boundary 
conditions to spin-dependent interfaces by taking differ- 
ent transmission ensembles for the two spin directions. 
These ensembles can differ in the total number of chan- 
nels and/or in the transmission distribution. Thus, we 
replace the spin-dependent current through the interface 



2e 
G A (e) 
4e 



(30) 



(l-/f( e ,0) + /^(-e,0)) , 



In general the spin-dependent quasi-particle conduc- 
tances Gf(e) and G^(e) entering the first term are of 
different magnitude and have different energy depen- 
dences. Similar as for the ferromagnetic wire we intro- 
duce the total conductance of the boundary Gs(e) = 
G-j-(e) + G|(e) and a dimensionless factor jB(epsilon) = 
(G-f(e) — G|(e))/GB(e), which we call polarization of 
the boundary conductance. Since the definitions ( p5| ) 
of quasi-particle and ( pfj| ) of Andreev conductance have 
been derived for a spin-degenerate interface these defini- 
tion are not valid anymore for spin-dependent interface 
scattering. It is, however, reasonable to assume that the 
energy dependence of all conductances is well approxi- 
mated by the same transmission ensemble, but different 
numbers of channels. We can motivate this choice by the 
fact that in the experiments that we have in mind the in- 
terfaces are strongly disordered regions, with a possible 
formation of an alloy layer extending over several mono- 
layers. In such contacts the number of channels is more 
or less controlled by the differences of the cross sections 
of the Fermi surface. But, on the other hand, the trans- 
mission ensemble and, hence, the energy dependence of 
the conductance is not expected to vary much in typical 
disordered contacts on the scale of the superconductor 
gap. 




FIG. 3. Resistance change of different types of contacts be- 
tween a normal metal reservoir and a superconducting reser- 
voir. The mixed contact (solid) varies from ballistic (q — 0) 
to tunnel (q = 1) from bottom to top. Intermediate values 
are (q = 0.25, 0.5, 0.75). The diffusive contact is shown by the 
short-dashed line and the dirty contact by a long-dashed line 
(see Table Q for a definition of these contacts). 
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We will in the following only take into account the 
differences in magnitude, but not in energy dependence. 
In the language of transmission distributions this means 
that the distributions are the same, but the number 
of channels differ. In this approximation the spin- 
polarization of the boundary conductance 7b is energy 
independent. The energy dependence of the Andreev 
conductance follows from the same transmission ensem- 
ble, but its magnitude will be reduced in comparison to 
the unpolarized case. It is important to notice that the 
boundary polarization and the polarization of the ferro- 
magnetic wire need not to have the same sign, since they 
are parametrically independent. The possibility of this-is 
demonstrated by microscopic numerical calculations B2l 



IV. RESULTS AND DISCUSSIONS 

We solve the kinetic equations presented in Sec. |l| in 
the three cases: 

A. purely elastic scattering 

B. inelastic scattering in linear response 

C. inelastic scattering in nonlinear response 

The boundary condition on the superconducting side of 
the wire has been derived in Sec. III. The boundary con- 
ditions at the ferromagnetic reservoir are 



/ T (e,-L) = / i (e,-L) = /(e;eF,T) 



(31) 



where f(e;eV,T) = (exp((e - eV)/k B T) + l)" 1 is the 
Fermi-Dirac equilibrium distribution at a constant volt- 
age V and temperature T. 

In the case of inelastic scattering Eq. ([H]) also implies 
that the electron temperature equals the lattice temper- 
ature in the ferromagnetic reservoir T e \(x = —L) = T. 
The other boundary condition for the electron temper- 
ature comes from the conservation of energy current in 
the ferromagnet and into the superconductor. 

As a reference we calculate the resistance of the system 
in the normal state 



-RpN — R% 



R 



BN 



(7 - 7b) 2 

1 + i? s f / i?BN 



(32) 



The third term is due to the spin-accumulation in the fer- 
romagnetic wire determined by the 'spin-flip resistance' 
R si = 1/G sf = R F l sf /Lta,nh(l B {/L)(l - 7 2 ). In the limit 
of a weak ferromagnet 7 2 <C 1 and a short spin-flip re- 
laxation length Z s f <C L the spin-flip resistance reduces 
to R s i ~ Rplsf/L, i.e. the resistance of a piece of the 
ferromagnetic wire of length l s {. We see that the excess 
resistance due to the spin-accumulation increases with in- 
creasing asymmetry between the polarization of the bulk 
conductivity and the polarization of the interface conduc- 
tance. The expression (|32|) will be used in the following 



to calculate resistance changes below the transition to 
the superconducting state: 



AR FS (T, V) = Rfs(T, V) - Rfn 
The differential resistance is defined by 



(33) 



-Rfs 



In the linear response regime we will omit the arguments 
of the differential resistance i?ps = Rfs(T, V — > 0). 

In the following analysis it will be useful to define the 
following temperature dependent average 



df(e;0,T) 
de 



de 



(35) 



This average occures, e.g., in the temperature dependent 
conductances of a contact between an incoherent metal 
and a superconductor (^9|). 



A. Elastic scattering 

When the scattering in the wire is elastic, the general 
solution of (m satisfying (pL) may be written as 



Me, x) = (/(e; 0, T) - f(e; eV, T)) (l + |) 



G F 

+f{e;eV,T) 



(36) 



The spatially independent spectral conductance Gfs(e) 
determines the current through the structures and re- 
mains to be found. The solution of the second kinetic 
equation (^|) satisfying the boundary condition ( |3l| ) is 



f sp (e,x) = 2asinh 



L + x 
I, 



(37) 



where the parameter a should be found from the continu- 
ity of the spin currents into the superconductor (|3l|) and 
ferromagnetic wire (0). We find the electrical current 



I(T, V) 



1 

Ye 



deG FS (e) x 



[l-/(e;eV > T)-/(-e;eV 1 T)] . 



(38) 



This expression shows that the spectral conductance de- 
termines the transport in each energy slice depending on 
the difference in occupation of states at this energy in 
the reservoirs. This form is analogous to the classical 
definition of a conductance as the proportionality factor 
between current and voltage difference. 
The spectral conductance is given by 



1 



1 



1 



G FS (e) G F G QP (e) + G A (e) 



(39) 



' B G Q p(e)+G A (c) 



G s t + Gqp(c) ( 1 — 7b 



»(<0 



:Gq P (e)+G 4 (e), 
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In the general case the full expression has to be used to 
calculate the resistance change in the superconducting 
state. 

When the ferromagnetic wire dominates the resistance 
of the whole structure a simplified expression for the lin- 
ear resistance change may be obtained. We first limit the 
discussion to the case of a weak ferromagnet and vanish- 
ing boundary polarization to obtain 



Ai? FS = 



1 



G Q p(e) + G A (e) 
1 



— -Rbn 



(40) 



Gsf + Gqp (e) / G s ( + G 



'BN 



We see that the resistance change consists of two con- 
tributions. The first is the resistance change due to the 
change of the boundary resistance, which would also be 
present in the absence of spin polarization. Note, how- 
ever, that this term can be qualitatively different from 
the case of a normal metal wire in contact to a supercon- 
ductor since in this case the proximity effect would not be 
negligible. The second term accounts for the difference in 
spin accumulation between normal and superconducting 
state. 

First we discuss the influence of spin accumulation 
on the FS-resistance for a spin-degenerate interface. 
In Fig. ^ resistance changes for two types of contacts 
are shown for different polarizations of the ferromag- 
net. Solid curves are for a relatively good contact 
(q = 0.75) and dashed curves for a less transparent con- 
tact (q — 0.25). In this plot the total resistance of the 
system is dominated by the resistance of the ferromag- 
netic wire i?F = 100i?BN and the spin-relaxation length 
is l s t = 0.03L, resulting in a spin accumulation resistance 
R s f w 3i?BN- Accordingly, the resistance change is nor- 
malized to i?BN to show the relevant scale of the effect 
produced by the superconducting transition. For both 
contacts spin accumulation (increasing from the bottom 
to the top curves) leads to an enhancement of the resis- 
tance. Specifically the low temperature resistance is well 
accounted for by Eq. (Eol) in the limit T — > 0: 



AR FS (T = 0) 



1 



G A (0) 



- i?BN + r 



R 



st 



1 + RBN/Rsf 



(41) 



The second term of this equation shows that the spin- 
accumulation always enhances the resistance, maximally 
by an amount 7 2 i? s f- The enhancement for the q = 0.25- 
contact has a uniform temperature dependence and does 
not change qualitatively. This is different for the q = 0.75 
contact. Here the resistance decreases monotonically in 
the unpolarized case as a result of the Andreev enhanced 
conductance. A small polarization 7 w 0.2 — 0.4 results 
in a nonmonotonic temperature dependence, i.e. an in- 
crease of resistance at lower temperatures. This can lead 
(for specific parameters) to a reentrant behavior of the 
resistance change, even overshooting the normal state 
value for larger spin-accumulation. At even higher spin 



polarizations 7 2 > Rbn/Rs! the Andreev contribution is 
completely masked and the resistance increases monoton- 
ically. This behavior resembles that of a less transparent 
contact if the absolute scale is properly chosen. 



< 




FIG. 4. Temperature dependence of the resistance change 
of the F-wire attached to a superconducting reservoir. Re- 
sults are presented for two mixed contacts with q — 0.5 (solid 
lines) and q = 0.25 (dashed lines). The conductivity polariza- 
tion 7 of the wire is changed from to 0.6 in steps of 0.2 from 
the bottommost curve to the topmost curve for both con- 
tacts. Other parameters are l s f = 0.03L, and Rf = 1007?bn- 
Clearly spin accumulation leads to an enhanced resistance in 
both cases. The resistance of the (/=0.25-contact is more or 
less uniformly increased. The effect of spin accumulation is 
much more dramatic for the g=0.75-contact. The monotonic 
resistance decreases in the unpolarized case is first turned 
into an reentrant behavior for small polarization overshoot- 
ing the normal state resistance slightly at low temperatures. 
Increasing 7 further leads to an increased resistance for all 
temperature. Note that this behavior resembles that of a 
q=0.25-contact, if properly rescaled. 



Let us now discuss the effect of the interface polar- 
ization 7b on the resistance change. In Fig. ^| the tem- 
perature dependent resistance of a q — 0.75 contact is 
shown for different interface polarizations. Other pa- 
rameters are Rf = 100i?BN, hi = 0.03L, and 7 = 0.3. 
The interface polarization 7b changes from the symmet- 
ric value +0.5 to the antisymmetric value —0.5, as indi- 
cated in the plot. The reduction of the Andreev conduc- 
tance by the spin-dependent interface resistance is taken 
into account by a phcnomcnological renormalization fac- 
tor (1 — 7b). To gain some insight it is useful to look at 
the low temperature limit of the resistance change in the 
limit J?bn < -Rsf- From Eq. (p9|) it follows that 



Ai? FS (T = 0) = 



G A (0) 



- i?BN + i? s f(4 7 7B - 7b) • ( 42 ) 



The spin-dependent contribution depends on the rela- 
tive sign of the two polarizations 7 and 7b and can also be 
negative (if 4773 < 7b)- This effect is seen from the lower 
two curves in Fig. ^ with an antisymmetric interface po- 
larization. An increasing interface polarization leads to 
a lowering of the resistance change, despite the increase 
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of the resistance due to the renormalization of the An- 
dreev conductance. It is worthwhile noting that for the 
largest negative interface polarization shown (73 = —0.5) 
the total resistance drop is larger than the resistance 
drop which would result from the pure Andreev reflec- 
tion in the absence of spin polarization of the interface 
and the F-wire. This apparent contradiction to the in- 
tuition that any spin-accumulation should decrease the 
Andreev-caused resistance drop stems from the fact that 
we plot the resistance change below the superconduct- 
ing transition. The contradiction is resolved by noting 
that the total resistance Rps(T) = i?FN + ARps(T) is 
always higher than for the unpolarized case. However, in 
a real experiment (with fixed polarizations) the Andreev 
conductance in the absence of a polarization can not be 
measured separately. It may therefore appear that the 
measured resistance drop is larger than one would expect 
from a simple estimate of the reduction of the interface 
resistance due to Andreev reflection. 



An analytical expression for the total conductance of 
the system can be found in the linear response regime. 
In this regime the effects of electron heating vanish since 
they will only contribute to the current in higher or- 
ders of the source-drain bias. The coupled equations 
for the spin-dependent chemical potential distributions 
and the electron temperature are simplified by letting 
T e i(x) -> T. By solving (M) and (O) together with 
the boundary condition ([23|) and ([31[) we find the lin- 
ear response resistance. Assuming a weak ferromagnet, 
7 2 <C 1, and a small interface resistance compared to 
the resistance of the ferromagnetic wire i?BN *C Rf the 
resistance change can be written as 



Ai?F S (T) 



1 



<G Q p(e) + G A (e)> 



R 



BN 



(43) 



[G s f + (Gqp(e)) G s[ + G 



BN 




-0.2 



-0.4 : 



0.5 1 

T/T 

c 

FIG. 5. Effect of the relative polarizations on the resistance 
change. The contact is a mixed contact with q = 0.75, other 
parameters are 7 = 0.3, Z a f = 0.03L, and Rf = 100-Rbn- 
The polarization of the boundary conductance 7b is varied 
between a symmetric configuration 7b = 0.5 and an anti- 
symmetric configuration 7s = —0.5. At the same time the 
Andreev-conductance is rescaled by a factor 1 — 73 to ac- 
count for the smaller number and transmission of Andreev 
channels. For large antisymmetric polarization the resistance 
decrease exceeds the decrease of the corresponding normal 
metal-superconductor contact. 



B. Inelastic scattering - linear response 

We will now proceed to study the case of inelastic scat- 
tering in the ferromagnetic wire. It is assumed that the 
current in the ferromagnet is weakly polarized, 7 -C 1. In 
order to simplify the discussions we disregard the possible 
asymmetry in the interface transparency in the following 
discussions and set 7b = 0. An extension is straightfor- 
ward. 



The first two terms in (|43|) are due to the effective in- 
terface resistance between the ferromagnet and the su- 
perconductor and the third term is due to the spin- 
accumulation. The latter term vanishes when 7 — > or 
l s [/L — ► 0. This equation has to compared with Eq. ([33] ) 
for the case of purely elastic scattering. Only the quasi- 
particle conductance enters the spin-accumulation contri- 
bution since spins cannot be injected into the supercon- 
ductor by means of the Andreev process. The tempera- 
ture averaged conductances directly determine the tem- 
perature dependence of the total resistance in the case of 
dominant inelastic scattering processes. The quasi parti- 
cle conductance vanishes at zero temperature since then 
no spin-current can propagate into the superconductor. 
At zero temperature (Gqp(c) + Ga(c)) = Ga(0) and the 
resistance the FS system in the case of inelastic scatter- 
i ng ( j43| ) equals the result in the case of elastic scattering 

The results with inelastic scattering in general differ 
from those with purely elastic scattering when the tem- 
perature is non-zero or when the current is measured 
in the non-linear source-drain response regime. The re- 
markable difference between Eq. ( |33| ) and Eq. ([l3|) is the 
way the thermal averaging is carried out. E.g. in the first 
term we have to average the inverse contact conductance 
in the case of elastic scattering, whereas we first have 
to average the conductance and than invert the result in 
the case of inelastic scattering. A similar consideration 
holds for the spin-accumulation term. The origin of this 
difference can be understood in the following way: we 
may visualize our wire (or any system) as mapped onto 
an electric circuit which contains energy-dependent con- 
ductors. In the case of purely elastic scattering we first 
have to calculate the total conductance of the system 
for each energy. The current is then found by averaging 
this spectral conductance with the difference of distri- 
bution functions of the adjacent reservoirs. This proce- 
dure yields Eq. (B3) for the change in the resistance. In 
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contrast, inelastic scattering equilibrates the local distri- 
bution of electrons in a way that the chemical poten- 
tial is equal to the potential found from solving the cir- 
cuit problem of the corresponding electric circuit. Thus, 
Eq. (^) follows the Kirchhoff's laws for our system. As 
we will demonstrate below this difference can have sig- 
nificant consequences for the temperature dependence of 
the resistance. 

Let us now illustrate the temperature dependence 
of the linear response conductance in the case of a 
metallic contact with an interface conductance much 
larger than the conductance of the ferromagnetic wire. 
The total conductance at sufficiently low tempera- 
tures is then (G QP (e) + G A (e)) w G A (0) = 2G BN 
and the quasiparticle conductance is (Gqp(c)) ~ 
(87rfc B T/A) 1 / 2 exp(-A/fc B T). The temperature must 
then be so low that 



/c R T< 



A 



ln(i?BN/i?sf) 



(44) 



in order to prevent thermally assisted spin-current into 
the superconductor. 
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FIG. 6. The ratio of the resistance Rfs (T) to the interface 
resistance -Rbn as a function of the reservoir temperature T. 
The ferromagnetic wire and is parameterized by the polariza- 
tion 7 = 0.3, the spin-flip diffusion length l s { — 0.2L and the 
interface resistance is Rbn ~ 0.05Rf. The spin-flip resistance 
is Rsf = 0.2R F . 

We show in Fig. |^ the ratio of the linear response re- 
sistance change A-Rps to the interface resistance Rbn 
as a function of the temperature T for a metallic in- 
terface with Rbn = 0.05-Rf, polarization 7 = 0.3 and 
spin- flip diffusion length l s f/L = 0.2. For these param- 
eters we have a 'spin-flip' resistance corresponding to 
R s f = 0.2-Rp. The change in resistance below the super- 
conducting transition temperature is due to a competi- 
tion between the excess resistance caused by the spin-flip 
relaxation and the reduced interface resistance caused by 
Andreev reflection. At T = we find from the approx- 
imate result ( j43| ) that i?Fs — -Rfn = 0.5-Rbn — 7 2 -R s f = 
-0.14.Rfn rougly corresponding to the numerical value 



which has been obtained without making the approxi- 
mation 7 2 <C 1 and i?p 3> Rbn ■ Using the condition 
( pd] ) we find that the spin-accumulation is strongly re- 
duced around T/T c = 0.7 and consequently the resistance 
of the system decreases before increasing again around 
T/T c = 1 where the boundary resistance is increased. 
This explains the non-monotonic behavior of the linear 
response resistance as a function of the temperature. 



C. Inelastic scattering - nonlinear response 

At a finite bias voltage the electron heating effects have 
to be taken into account and the coupled equations for 
the electron temperature, the spin-dependent chemical 
potentials ([l0|), (pd]) , ([I4 ), and ( |l4|) and the boundary 
conditions ( p3[ ) and (|3l|) have to be solved numerically. 

First let us discuss the transport properties when the 
electron-phonon interaction is weak so that we have per- 
fect conservation of energy current and the left hand side 
of ( |l4] ) can be set to zero. From the discussions in the 
previous section we understand that there will be a reduc- 
tion in the excess resistance due to the spin-accumulation 
when the electron temperature on the ferromagnetic side 
reaches condition (Q) so that there is a significant spin- 
current entering the superconductor. Roughly speaking, 
the electron temperature on the ferromagnetic side is pro- 
portional to the applied source drain bias. Thus, as a 
crude approximation, we expect that the excess resis- 
tance due to the spin-accumulation is lowered when 



iV < 



A 



In (i?BN/i?sf) 



(45) 



We show in Fig. || the resistance change Rfs(T = 0, V) — 
Rbn d33| ) normalized by the interface resistance Rbn as 
a function of the bias voltage V. As before, the interface 
resistance is -Rbn = 0.05i?F, the polarization 7 = 0.3 
and the spin-flip diffusion length l s f/L = 0.2. For these 
parameters we have a 'spin-flip' resistance corresponding 
to R s f = 0.2i?p- The change in resistance below the su- 
perconducting gap is due to a competition between the 
excess resistance caused by the spin-flip relaxation and 
the reduced interface caused by the Andreev reflection. 
A dip in the resistance is seen around V = 0.7 A which is 
correctly described by ([45|). Below this bias voltage the 
resistance is caused by the competetion between spin- 
accumulation which enhances the resistance and the ef- 
fective interface resistance which reduces the resistance. 
At higher voltages the resistance is only caused by the 
effective interface resistance and the reduction of the re- 
sistance change Rfs(T — 0, V) — Rbn as a function of 
the bias voltage is small. 

In the limit of strong electron-phonon interaction the 
electron temperature equals the lattice temperature. The 
spin-current into the superconductor is then not en- 
hanced due to thermal activation and consequently the 
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spin-accumulation on the ferromagnetic side is only re- 
duced when the potential on the ferromagnetic side of the 
interface is higher than the superconducting gap. This 
occurs when 



eV = All 



Rf 
Rbn 



(46) 



and thus at a potential that is much larger than the 
superconducting gap, in contrast to the case of weak 
electron-phonon interaction. 
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FIG. 7. The ratio of the resistance Rfs(T = 0, V) to the 
interface resistance -Rbn as a function of the bias voltage V . 
The ferromagnetic wire and is parameterized by the polariza- 
tion 7 = 0.3, the spin-flip diffusion length l s { — 0.2L and the 
interface resistance is Rbn = 0.05i?F. The spin-flip resistance 
is R si = 0.27? F . 

In the intermediate regime the electron-phonon in- 
teraction should be included. In order to illustrate 
the main physics we consider the case of a weak po- 
larization and set 7 = and consequently there are 
no effects due to spin-accumulation and the resistance 
change of the wire is only due to the change of the ef- 
fective boundary resistance. The chemical potential in 
the ferromagnetic wire is thus spin-independent. Fur- 
thermore we consider the case that the lattice temper- 
ature is zero, T = so that the electron temperature 
arises solely due to electron heating. Solving the diffu- 
sion equation (|l^) on the ferromagnetic side of the inter- 
face gives n{x) = —eVx/L + /i(0) [1 + x/L], where eV is 
the applied bias and fj,(0) is the potential drop across 
the ferromagnet-superconductor interface. The super- 
conducting energy gap A presents a natural energy scale 
for the problem. We will characterize the strength of 
electron-photon energy exchange by a dimensionless con- 
stant k = AL(e 2 A 3 /GF, C is defined by the relation (15). 
The energy diffusion equation then simplifies to 

2 

y (L T d x f (fc B T cl /A) 2 = (fc B T cl /A) 5 

-((eV--/i(0))/A) 2 /^ (47) 
where we introduce a typical length-scale for the energy 



exchange Lt = L/^/k. If n <C 1 L <C Lt and the ex- 
change is not effective. For longer wires, k becomes big- 
ger than unity. In this case, the electron temperature de- 
velops a constant plateau in the ferromagnetic wire and 
only changes rapidly within the length-scale Lt near the 
end-points x = —L and x = 0. It follows from (47) that 
in this case the temperature in the middle of the ferro- 
magnetic wire becomes 



k n T c 



1 _ 1/5 / eV- n(0) 



2/5 



(48) 



We will now present numerical results of the temper- 
ature profile in the ferromagnetic wire and the result- 
ing resistance change using 7 — 0, a metallic interface 
Rbn/Rf = 0.05 for various values of the electron-phonon 
coupling interaction. We show in Fig. || the spatially de- 
pendent electron temperature in the ferromagnetic wire 
for k = 10 6 at a bias voltage eV = 40 A (upper curve), 
eV — 20 A (mid curve) and eV — 10 A (lower curve). The 
electron temperature in the middle of the wire follows 
from ( [l8"| ) . There are rapid changes of the electron tem- 
perature close to the ferromagnetic and superconducting 
reservoirs and the temperature in the middle of the wire 
is lower than the electron temperature close to the super- 
conductor. The latter temperature is important for the 
effective interface resistance. 
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FIG. 8. The spatially dependent electron temperature in 
the ferromagnetic wire. The upper curve is for a bias voltage 
V — 40A, the mid curve for a bias voltage V = 20A and the 
lower curve for a bias voltage V = 10A. The metallic inter- 
face resistance is Rbn = 0.05-Rf, the reservoir temperature 
T = and the electron-phonon interaction strength k = 10 6 . 

We show in Fig. || the resistance change as a function 
of the bias voltage. The different solid lines show the cur- 
rent for different ratios of the electron-phonon coupling 
starting from no electron-phonon interaction (a) k = 
going through intermediate electron-phonon interaction 
(b) k = 10 2 , (c) k = 10 6 , (d) k = 10 8 to strong electron- 
phonon interaction (f) k = 00 when the electron tem- 
perature equals the lattice temperature, e.g. when there 
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is no energy transfer between the electron and phonon 
system. The cross-over bias voltage for the excess resis- 
tance is sensitive to the strength of the electron-phonon 
interaction and occurs from around A (a) to around 40A 
(f) (according to (fi"6|)). The dependence on the electron- 
phonon interaction parameter k is rather weak as can be 
understood from (fi"8|). The local electron temperature in 
the middle of the ferromagnetic wire is proportional to 
k -1 / 5 and thus only has a very weak dependence on n. 
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FIG. 9. Resistance change (Rfs(T = 0, V) — Rfn) / Rfn as 
a function of the bias voltage V. The ferromagnetic wire 
is described by the polarization 7 = 0, and the interface 
resistance is Rbn = 0.05Rf. Curve (a) corresponds to no 
electron-phonon interaction k — 0. Curve (b), (c), (d), and 
(e) correspond to intermediate electron-phonon interaction, 
k — 10 2 , K = 10 4 , K = 10 6 and k — 10 8 , respectively. 
Curve (f) is for the case of strong electron-phonon interac- 
tion K — OO. 



V. DISCUSSION OF EXPERIMENTS 

In this section we discuss the connection of our results 
with experiments of Petrashov et alB and Giroud et al.a. 
It will turn out that most of the experimental results can 
be understood on the basis of our calculations. Both ex- 
perimental arrangements we will discuss below contain 
F-S junctions where the superconductor and the ferro- 
magnet overlap in a certain region. The current redis- 
tribution in these junctions will play an important role 
in the following. Let us therefore introduce parameters 
characterizing these juntions: the resistance of the in- 
terface is called -Rbn in accordance with our previous 
consideration. Additionally i?sj will be the resistance of 
the superconducting part of the overlap junction in the 
normal state and Rvj the resistance of the ferromagnetic 
part of the overlap junction. 

In the experiment by Giroud et ala a non-monotonic 
behaviour of the resistance below the superconducting 
critical temperature was observed. The sample consisted 



of a ferromagnetic wire, the resistance of which was mea- 
sured in a 4-point arrangement. At some point a super- 
conducting strip was on top of the wire. In a second 
sample two such strips were present and the resulting re- 
sistance change was twice as big as in the case of one 
strip. Since our formulation is based on a single interface 
and no coherent coupling between the two superconduct- 
ing strips was found experimentally, we concentrate here 
on the sample with one strip. The resistance change in 
the two-strip sample is then simply twice that for the 
single strip sample. The experimental arrangement is 
such that in the region of the strip the current is redis- 
tributed among the ferromagnet and the superconduc- 
tor. In Appendix ^ we introduce a simple quasi-one- 
dimensional model to calculate the effective resistance, 
these results being used for comparison with experiment. 
The resistance of the superconducting Al-strip is QAil, 
the resistance of the ferromagnetic part below the strip 
is 10O, and the resistance of the interface is estimated to 
be O.lfi. Since the measured resistance change of the F- 
wire shows no signature of the vanishing of the resistance 
of the superconducting part, we believe that the real in- 
terface resistance is higher than estimated in Ref. ^, in 
particular higher than i?s j ■ This yields a total resistance 
of i? e ff = 2(i?pji?BN) 1 ^ 2 , which is approximately of the 
order of a few fl. A resistance change of the interface re- 
sistance Ai?B will than lead to an change of the effective 
resistance AR c s = (Rpj / Rbn) 1 / 2 ARb, which in the case 
i?Fj > ^?bn is large than the resistance change of the in- 
terface resistance itself. For the experimental values we 
have (Rfj/Rbn) 1 ^ 2 ~ 10 and thus a resistance change of 
« 0.2f2, as observed in the experiment may result from a 
change of the interface resistance JXbn ~ O.lfi by 20%. 

The results of Petrashov et alu are more intriguing, 
since the magnitude of the measured resistance drop in 
some of the samples seems to be far too large to be ex- 
plained without a 'long-range' proximity effect in the 
ferromagnet. We will concentrate here on three of the 
four samples discussed in Ref. ||. In these samples the 
transport through a long ferromagnetic wire with one fer- 
romagnetic and one superconducting contact is studied, 
this is in contrast to experiments of Ref. |^. The geometry 
is such that the superconducting contact overlaps the fer- 
romagnetic wire at one end and the current has to pass 
through a tiny piece of the superconductor. The three 
samples differ in the interface resistance. Two samples 
with a low interface resistance show large drops of the 
resistance of the order of 8f2 respectively 16S1 below the 
superconducting critical temperature. The third sample 
has a higher interface resistance (-Rbn = 41f2) and shows 
a small resistance increase of the order of 1.50. 

This agrees qualitatively with the results of our model. 
Indeed, the bigger resistance of the boundary usually 
means a formation of a ticker tunnel barrier such that the 
transmission eigenvalues are shifted towards zero. Our 
model does predict a resistance decrease for a fairly trans- 
parent interface and changes to an increase for a more 
tunnel-like interface. This is shown e.g. in Fig. El. How- 
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ever, quantitatively one would expect that the resistance 
changes below the superconducting transition tempera- 
ture are always of the order of the boundary resistance 
itself. 

This is obviously not the case in the experiment with 
-Rbn = 410, where the measured resistance change is 
about 40 times smaller. The first idea is that the re- 
sistance drop of the samples with better interface may 
possibly be accounted for by combining the effect of 
current redistribution and the apparent enhancement of 
Andreev reflection discussed in Sec.jfy]. Again we cal- 
culate the effective interface resistance in a quasi-one- 
dimensional model (see Appendix ^|). In the limit of 
a small interface resistance i? e ff2 = (RfjRbn) 1 ^ 2 ■ A 
change of the interface resistance again results in an ap- 
parently larger change of the effective resistance AR e ff2 = 
(Rfj/RbnY^SRb/^. We may speculate that the large 
resistance drop observed in the experiment by Petrashov 
et al. can possibly be explained by this effect together 
with the observation made in IV that a spin dependent 



interface may cause another apparent enhancement of the 
Andreev- reflection. 
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FIG. 10. Comparison of experiment of Giroud et a/.i and 
theoretical calculations. The two experimental data sets are 
for the sample with one superconducting strip (squares) and 
for the sample with two superconducting strips (circles). The 
theoretical curves are obtained from Eq. ( p9| ) with the fol- 
lowing parameters: i?p = 1007?bn, spin relaxation length 
Isi = 0.012, conductivity polarization 7 = 0.3 and a double 
interface contact. All curves are normalized to the respective 
maximal value. 



There may be a more radical explanation for a small 
relative resistance change. In fact, the morphology of 
the metal-ferromagnet interfaces has not been yet suf- 
ficiently studied. The. actual structure of the interface 
may be complicated.E2l To illustrate how this can affect 
the results let us consider a simplistic model of a double 
interface. We speculate that a thin layer of magnetic alloy 
separates the ferromagnet and the superconductor. The 
boundary scattering then occurs in two stages: at the 
"inner" interface between the ferromagnet and the alloy 



and at the "outer" interface between the superconduc- 
tor and the alloy. Since the proximity effect is quenched 
in magnets, the resistance of the "inner" interface is not 
affected by the superconducting transition whereas the 
resistance of the "outer" interface acquires a change de- 
scribed above. This leads to a smaller relative resistance 
change. 

In Fig. [l^ we show a comparison between the exper- 
imental results of Giroud et al. and our calculation for 
a contact with a transmission eigenvalue distribution of 
the model just described (see Appendix |b|). The maxi- 
mal relative resistance change of the contact is ss 16% in 
our calculation. According to above considerations this 
may be enhanced to a measured resistance change of the 
order of w (i? F j/^BN) 1/2 16%(5i?B ~ 160%<5i? B in agree- 
ment with measurements. 



VI. CONCLUSIONS 

To summarize, we have calculated the resistance of 
a diffusive ferromagnetic wire in contact with a super- 
conducting reservoir in the linear and non-linear regime 
with purely elastic and inelastic scattering. It has been 
demonstrated that most of the recent experimental re- 
sults can be understood in the absence of a supercon- 
ducting proximity effect in the ferromagnet. Spin accu- 
mulation leads to an enhanced resistance below the su- 
perconducting transition temperature whereas Andreev 
reflection can lead to a decreased resistance below the 
superconducting transition temperature. The competi- 
tion between these two mechanism determines the sign 
of the resistance change. The magnitude of the resis- 
tance change is of the order of the interface resistance or 
the spin-relaxation resistance. Electron heating can dra- 
matically modify the nonlinear response resistance and 
change the bias dependence by orders of a magnitude. 
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APPENDIX A: CURRENT REDISTRIBUTION IN 
AN OVERLAP JUNCTION 



Here we introduce a quasi-one-dimensional model to 
account for the redistribution of the current under an 
overlap junction as used in the experiments. The two 
geometries we have in mind are depicted in Fig. Ell 
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FIG. 11. Schematic picture of current redistribution in an 
overlap junction. The geometry a) corresponds to the exper- 
iment by Giroud et alB. If the resistance of the ferromagnet 
below the overlap region is the highest, a considerable part 
of the current flows through the interface and the supercon- 
ductor. The gesmetry in b) corresponds to the geometry of 
Petrashov et ala. Here the current is forced to leave the con- 
tact region through the superconductor. 



To estimate the measured resistance we use the follow- 
ing quasi-one-dimensional model for the current redistri- 
bution in the overlap region of length d. The currents 
If(x) in F and Is(x) in S in direction of the ferromag- 
netic wire follow from Ohm's law 



If(x) 



d dUp(x) 
i?F,i dx 



d dU s {x) 
Rsj dx 



(Al) 



where -Rfj(sj) is the resistance of the ferromagnetic part 
under (superconducting part above) the contact and 
Uprg) the respective voltage. Current conservation dic- 
tates 



dI F (x) dl s (x) U s (x)-U F (x) 



dx 



dx 



R B Nd 



(A2) 



is the current per unit length through the contact resis- 
tance is Rbn- Boundary conditions are obviously that 
the total voltage drop is equal to V and no current leaves 
the system through the boundary to vacuum. 

Solving these equations for the geometry of Giroud et 
we find the effective resistance of this part to be 



Re 




(A3) 



-Rsj + Rfj 



tanh 



'-Rfj + #SJ 
4i?BN 



Of specific interest if the case that -Rfj ^> Rbn 3> -Rsj 
in which case we obtain i? e ff w 2(i?pj-RBN) 1 ^ 2 - 

A calculation similar- to the previous for the geome- 
try of Petrashov et alU leads to an effective measured 
resistance of 



R 



V -RfjRbN 



■cS2 



tanh y/ Rfj / Rbn 



(A4) 



in the limit of vanishing resistance of the superconductor 
on top of the ferromagnet. The difference to the previous 
calculation is that here the current enters the junction 
through the ferromagnet, but has to leave the junction 
through the superconductor. 



APPENDIX B: TRANSMISSION EIGENVALUES 
OF A DOUBLE INTERFACE 

The distribution of transmission eigenvalues of a dou- 
ble interface as described in the text can be found with 
the technique described in Ref. |2^. For details we refer 
to these articles. We model the double interface by a 
ballistic contact and a tunneling barrier in series. The 
tunneling barrier of conductance Gt models the sharp 
drop in the potential due to the band structure mismatch, 
whereas the region close to that interface is treated as a 
collection of unit transmission channels with a total con- 
ductance Gqpc The distribution of transmissions can be 
found from the solution of 



/($) = G T sin($ 



2G 



QPC 



tan 



(Bl) 



The distribution of transmission eigenvalues p(T) is 
found by analytic continuation into the complex plane 



1 



1 



-.Re 



2i acosh 



1 



(B2) 



The dependence on the two separate conductances may 
be eliminated in favor of the total conductance of the 
contact Gbn = GtGqpc/(Gt + Gqpc) and the ratio of 
the two a = Gt/2Gqpc- The transmission eigenvalue 
distribution then only depends on a. It is plotted in 
Fig. |lj for several values of a. For small values of a the 
contact is dominated by the tunnel barrier resulting in a 
shift of the transmission eigenvalues to lower values and 
a gap above a certain T. Higher a shift the distribution 
to larger transmission eigenvalues and a gap opens up fol- 
low transmission eigenvalues. For a range 0.1 < a < 0.5 
the distribution restricted to a finite interval of trans- 
mission eigevalues. At even higher values of a the upper 
gap closes and the distribution becomes more and more 
peaked at T = 1. 




FIG. 12. Transmission eigenvalue distribution for the dou- 
ble interface. The ratio a — Gt/2Gqpc is varied between the 
different curves. 
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